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Abstract 

A dynamical cross-over regime is revealed when exposing a classical two-dimensional ordered Josephson junction (JJ) array 
to evanescent waves and tuning the incident microwave power. At the lowest possible temperature for these experiments, 1.1 
K, and at the lowest power setting, -55 dBm, evanescent waves are transmitted without loss and the resonance exhibits a 
quality factor of "4200. A second, smaller resonance, which evolves with increasing power from the main resonance, is also 
investigated. In contrast to the behavior of the main resonance, this second peak grows as the incident power is increased and 
does not maintain a fixed resonant frequency for temperatures less than the superconducting critical temperature of niobium. 
The tunability of both resonances is studied as a function of temperature and microwave power. Finally we speculate that this 
dynamical crossover regime is evidence of a transition between two states of phase coherence where at low microwave power 
the JJ arrays are phase locked and at high microwave power the JJ arrays are unlocked.TJ^K . 



Negative index of refraction metamaterials, materi- 
als that simultaneously have a negative permittivity 
and negative permabihty, bend electromagnetic waves in 
counterintuitive directions. [T] One promise of these ma- 
terials is their ability to transform decaying evanescent 
waves into amplifying waves and, by doing so, extend 
the resolution of optical microscopes to sub-wavelength 
range. ^ However, this opportunity [3] is jeopardized 
by dissipative losses which grow as the size of the con- 
stituents, the meta-atoms, are made smaller. 0] An en- 
couraging remedy to overcome losses is to directly em- 
bed gain media into metamaterials having fishnet geome- 
tries. Theoretical models suggest that as a result of loss 
compensation, electromagnetic signals can be amplified 
with 'high' quality factors, as 'high' as 800, for realistic 
values of gain. [5] While amplifying propagating electro- 
magnetic waves has successfully been demonstrated, al- 
beit with signals having a low quality factor, Q < 50, 
[S] a similar, tunable amplification of 77.on-propagating 
evanescent waves with high quality factors remains to 
be seen. And since losses and lack of tunability are the 
real show stoppers for future applications, achieving the 
remarkable properties known to metamaterials demands 
removing nearly all losses [7j while embedding nonlinear 
elements for tunability. 

Lossless materials which combine fishnet, or square 
loop, geometries with superconducting nonlinear gain 
media are known as two dimensional ordered Josephson 
junction arrays. Josephson junction (JJ) arrays have in- 
trinsic capacitances and inductances built into them - a 
key component in the design of negative index of refrac- 
tion metamaterials. As attractive tunable metamaterials 
their salient features are their low loss nature, nonlinear 
inductance, and miniaturized device elements. [5] All of 
this leads to materials that are well suited as metamate- 
rials but in practicality, JJ arrays, at least conventional 
JJ arrays, would not be of convenience since they need 



to be kept cold and radiated with microwaves. Never- 
theless, if this is done, then these are ideal materials for 
restoring decaying evanescent waves. 

In this Letter, we report low temperature measure- 
ments of evanescent wave amplification with quality fac- 
tors as high as 4200 when using unbiased superconduct- 
ing JJ arrays. Evanescent wave amplification, as de- 
scribed in the transmission experiments of Ref j£j , is the 
sharp increase in transmission through a waveguide with 
a metamaterial in the same frequency range as the cutoff 
regime of an empty waveguide. We also observe tun- 
ability and a surprising transition between two different 
sets of resonances with drastically different responses to 
temperature and incident power. Once the temperature 
exceeds the superconducting transition temperature, the 
two tunable resonances no longer prevail suggesting their 
relevance to the sample's superconducting properties. 

The samples are two dimensional JJ arrays that were 
commercially fabricated by the Sperry-Univac Corpora- 
tion, now Unisys, using the selective niobium annodiza- 
tion process [TU] and constructed on a silicon substrate 
covered with a 1 /im thick oxide layer. |11| These ar- 
rays, arranged in a square ring geometry are unique 
in that there is no ground plane. Each ring has four 
equally spaced Josephson junctions and is defined here 
as a "meta-atom" whose size is much smaller than the 
applied microwave wavelength, thus creating an effective 
medium. The length of one side is 29 microns, with a 
linewidth of 5 microns and a niobium thickness of 500 
nm. The junctions have an area of (2.5 microns)^ and 
are made of Nb-amphorous Si-Nb, with a-Si serving as 
a 10 nm thick insulating barrier. A portion of the JJ 
arrays is shown in Fig. la. The critical current of an 
individual junction is estimated to be 1.65 x 10~^ A. 
Transport measurements in a magnetic field on these ar- 
rays are reported elsewhere. [12] The superconducting 
transition temperature is T^, = 9.2 as measured by 



SQUID magnetometry. 

The original sample with one million (1000 x 1000) 
junctions is divided into several pieces with each piece 
containing only the arrays and not electrodes or pads for 
electrical contacts. The samples are placed in a Ag-plated 
Cu X-band rectangular waveguide and held in place with 
a piece of Rohacell, a styrofoam-like material with a di- 
electric constant between 1.06 - 1.11 and oriented such 
that the rf magnetic field is perpendicular to the surface 
of the sample. 

Attached to the rectangular waveguide are heaters, 
a thermometer and two antennas which are housed in 
a vacuum can inside a magnetically shielded cryostat. 
Semi-rigid coaxial cables [13] connect to the antennas 
and extend to a HP 8722 D vector network analyzer via 
phase-maintaining coaxial cables. |14) Standard trans- 
mission and reflection measurements are made by acqur- 
ing scattering parameters; the data is averaged over 16 
subsequent runs with 1601 data points taken for each 
run. No additional filters or amplifiers are added to the 
system. The cutoff frequency for an empty waveguide is 
6.56 GHz as shown in Fig. lb. 

We place several small, unbiased samples of different 
sizes in the waveguide. At room temperature, well above 
Tc of the J J arrays, the electromagnetic transmission has 
broad ripples and a broad resonance at 12 GHz as shown 
in the blue curve of Fig. Ic. Additionally, there is notice- 
able insertion loss of -7 dB and the signal from 5 GHz to 
the cutoff frequency is much broader than for the empty 
waveguide. By cooling to liquid helium temperatures, 
where the arrays are superconducting, several sharp res- 
onances appear at discrete frequencies both above and 
below cutoff as shown in the red curve of Fig. Ic. Also, 
the insertion loss goes to zero indicating that the lossy 
spectrum at room temperature is a consequence of the 
normal state JJ arrays and not the substrate. 

To explore these resonant features in greater detail and 
without the possibility of crosstalk, we place a single sam- 
ple of 0.5 cm x 0.5 cm dimensions with roughly 30,000 
meta-atoms at the center of the waveguide. At fixed tem- 
perature, T = 1.1 K, two distinctive transmission features 
appear below cutoff with increasing incident power from 
-55 dBm to - 5 dBm having drastically different behav- 
iors: the dominant resonant peak decreases in amplitude 
at fixed frequency then slightly shifts to lower frequen- 
cies over a range of ~ 0.4 MHz, until a second, smaller 
resonance develops and increases in amplitude with in- 
creasing frequency over a range of ~13 MHz as shown in 
Fig. 2a. Increasing the power also has the effect of broad- 
ening the dominant resonance, Q decreases from 4200 to 
340, while sharpening the smaller resonance, Q increases 
from 340 to 1240. The amplitude of the dominant res- 
onance is tunable over 30 dBm with a "net gain" of 20 
dB, and the smaller resonance is tunable over 10 dBm 
with a "net gain" of 8 dB. Here "net gain" refers to the 
increase in the signal's amplitude from negative dB to 
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FIG. 1: (color online) (a) Optical microscope image of a 
portion of the sample. Each square loop contains 4 Joseph- 
son junctions; the scale bar denotes 30 fj,m.{h) Plot of trans- 
mission vs frequency for an empty waveguide, (c) Plot of 
transmission vs frequency of the same waveguide as above 
but containing several unbiased samples at room (blue curve) 
and liquid helium (red curve) temperatures. Note: The value 
of -1-3 dB between 7-8 GHz for the liquid helium curve cor- 
responds to dB; this offset is due to the calibration being 
done at room temperature while the measurements are made 
at liquid helium temperatures. 



dB; this is not a real gain since the signal never exceeds 
dB. 

The observed transmission of evanescent waves from 
one end of the waveguide to the other occurs by tun- 
neling since evanescent waves cannot carry energy. The 
signature of 'transmission by tunneling' is equal, or near 
equal, amplitudes of forward and backward evanescent 
modes [S],[in]; this is shown in the reflection and trans- 
mission coefficients and plotted as a function of frequency 
for -55 dBm and -5 dBm in Fig. 2b. [16^ 

The resonant frequency, /, of a single JJ meta-atom is 
estimated from the geometric inductance, Lg, the Joseph- 
son inductance, Lj, and the capacitance, C, of the tun- 
nel junctions. This is expressed as f— ly — , "'" 



with Lg = [ill ( ~ i] where R is the radius of the 
square loop and r is the radius of the wire. [17J From the 
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FIG. 2: (color online) Measurements at T = 1.1 K. (a) Fre- 
quency dependence of array transmission with input power 
ranging from -55 dBm to -5 dBm, every -5 dBm. The arrows 
indicate increasing input power, (b) Magnitude of the trans- 
mission I S21 I (black circles) and reflection 1 Sn 1 (red circles) 
at -55 dBm and - 5 dBm. [16] 



dimensions of an J J meta-atom and the Josephson induc- 
tance which is defined as Lj — 2^7- with = ^ where 
Ic is the critical current, h is Planck's constant, and e 
is the charge of the electron, the estimated resonant fre- 
quency is ~ 50 GHz. This crude estimate is an order of 
magnitude higher than the observed resonant frequency 
suggesting that additional terms have been neglected, 
such as the kinetic inductance or possibly the estimated 
critical current is incorrect. In addition, this calculation 
is not accurate in describing the experimental situation 
since the resonant frequency depends on the location and 
orientation of the sample inside the waveguide. We com- 
pare this calculation to the Josephson plasma frequency. 




15 GHz, which is closer to the actual 



value. 

Given the discrepancies between predicted values and 
experimental observations, it is necessary to relate the 
observations to superconductivity and the Josephson ef- 
fect. For superconductivity, microwave measurements at 
different temperatures are made, with four different tem- 
peratures shown in color plots of Fig. 3. Perfect trans- 
mission, dB, at T =1.1 K and 6.0 K for -55 dBm is 
shown in green in Figure 3 a and b. 



Interestingly, as the temperature approaches T= 6.0 
K, higher and higher input powers are required for the 
emergence of the second peak; this behavior is clearly 
seen when comparing the transmission spectrum in Fig. 
3b to that in Fig. 3a. One possible explanation for this 
is that with increasing temperature, until T=6.0 K, more 
junctions become coherent with each other as larger num- 
bers of their individual linewidths overlap. Thermally in- 
duced coherence can be a result of the sample not being 
precisely centered in the waveguide or perfectly oriented 
with the rf magnetic field. The onset of coherence with 
temperature was previously mentioned in Ref. [TH] . 

This situation drastically changes for T > 6.0 K: the 
tunability of each resonance substantially decreases with 
increasing power. Indeed, since these are two dimensional 
JJ arrays, they are model systems for the Berezinskii- 
Kosterlitz-Thouless (BKT) phase transition. As previ- 
ously reported for these samples, the BKT transition 
temperature between the binding of vortex and anti- 
vortex pairs (T<Tbkt), and the unbinding of these pairs 
(T>Tbkt) occurs at Tbkt = 5.85 K. [H] Our data of a 
turning point at T= 6.0 K is consistent with this picture. 
At T = 8.0 K, the transmission of the smaller, higher 
frequency resonance saturates and is no longer tunable 
with increasing power as shown in Fig. 3c. However, the 
main resonance is still somewhat tunable at T=8.0 K, 
although transmission occurs at lower frequencies than 
the resonances at T ^ 6.0 K. At T= 10 K and beyond, 
there is only one broad, low amplitude resonance which 
is constant, albeit noisy, with incident power as shown 
in Fig. 3d; in aggregate, these temperature and incident 
power results are indicative of superconductivity. 

Central to understanding the nonlinearity and Joseph- 
son effect of these spectra requires explaining the tun- 
ability of the 6.314 GHz resonance for T ^ 6.0K. If the 
reduction in amplitude with increasing power is a result 
of dissipation, then as T increases, the amplitude should 
decrease. Keeping power fixed at -55 dBm and increas- 
ing T, the amplitude decreases. However, the decrease in 
amplitude at fixed power with temperature is not of the 
same ilk as the decrease in amplitude at fixed T with in- 
creasing power as seen in Fig. 2a and Fig. 3 . First, as T 
increases at fixed power, the resonant frequency shifts to 
lower values than for the case of increasing power at fixed 
T as shown in Fig. 4a. Secondly, and rather remarkably, 
tunability of the resonance occurs with very low incident 
power even though the niobium thickness is 500 nm. For 
even thinner films of niobium, ""200 nm, without JJs, sig- 
nificantly higher values of incident power are necessary, 
as high --1-20 dBm, to lower the amplitude by ~20 dB; 
for thin Nb films, this amplitude reduction, accompanied 
with a downward frequency shift, is due to dissipative 
losses. [1^ Lastly, and most surprisingly, is the decrease 
in loss, 1/Q, with increasing T for input powers between 
-40 dBm and -18 dBm for T < 6 K. This is shown in the 
shaded region of Fig. 4b and confirms that dissipation 
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FIG. 3: (color online). Microwave transmission (dB) as a 
function of frequency and input power for a.) T=l.l K, b.) 
T = 6.0 K, c.) T = 8.0 K, and d.) T =10.0 K. The color bar 
scale represents the range of transmission from zero transmis- 
sion (black) to perfect transmission (green). Note: The value 
of +4 dB corresponds to dB due to the calibration being 
done at room temperature while the measurements are made 
at liquid helium temperatures. 



is not driving the decrease in amplitude with increasing 
input power of the main resonance. 

If the tunability of the main resonance is not due to 
dissipation, then another possibility is that the JJs are 
operating as parametric amplifiers owing to their non- 
linear inductance. The similarity between these results 
and those due to parametric amplification might not be 
coincidental since the signal from the network analyser 
is not a pure tone, but rather contains sidebands. It is 
by these sidebands that we are likely seeing " gain" . We 
measured the output signal at 6.4 GHz from the network 
analyzer; there are sidebands spaced at -I-/- 20 kHz and 
-V j- 50kHz as well as higher order harmonics. In addition 
to a central frequency, sidebands from separate sources 
are standard for amplifying signals with JJ arrays; how- 
ever, in these experiments, the amplitude of the signal is 
always greater than dB.pP] - 

The losses, 1/Q, increase with increasing T for high 
power resonances, greater than -18 dBm, indicating in- 
creasing dissipation with temperature as shown in Figure 
4 b. Additionally, since these resonances shift to higher 
frequencies with incident power as shown in Figs. 2 a 
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FIG. 4: (color online), (a) Resonant frequency vs input power 
for different temperatures, (b) (1/Quality Factor) vs input 
power for different temperatures. The shaded gray region, - 
40 dBm to -18 dBm, indicates a regime where 1/Q decreases 
with increasing temperature from 1.1 K to 6.0 K. (c) 8^2 and 
Sfi vs input power for T=l.l K and T = 10 K. 



& 4 a this suggests that these resonances are not due 
to quasiparticles, but likely free vortices. When compar- 
ing 1/Q to the transmission and reflection coefficients we 
observe that as 1/Q increases or decreases, the reflection 
coefficients follow suit as shown for T = 1.1 K in the 
red curve in Fig. 4c. Above the superconducting tran- 
sition temperature, T = 10 K, the signal is all reflection 
as shown in the open circles of the blue curve of Fig. 4c. 
So, why then, does one peak evolve into a second peak? 
Why are there two peaks? What is the nature of the sec- 
ond peak? These are left as open questions and require 
further investigation. 

However, and without a corroborative theory to back 
up our findings, we naively speculate that the difference 
between these two sets of resonances is a measure of how 
phase sensitive the JJ arrays are to incident power: we 
speculate that the arrays are phase insensitive at low 
power and thus are phase locked for the dominant reso- 
nances, while for the second set of resonances, at higher 
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incident powers, the arrays are phase sensitive, and thus 
unlocked. This results in a constant resonant frequency 
for the phase locked arrays, but a shifting resonant fre- 
quency for the unlocked arrays with incident power. 

We also base our thinking on eliminating what these 
resonances are not: first, we rule out generation of quasi- 
particles with increasing incident power for T ^ 6.0 K be- 
cause in our work increasing power shifts the resonances 
to higher frequencies not lower frequencies as would oth- 
erwise be the case. Second, we rule out extraordinary 
transmission [23] since the resonant frequency and sam- 
ple location in the waveguide are coupled to each other. 

In conclusion, tunable amplification of evanescent 
waves has been demonstrated using unbiased, fishnet 
Josephson junction metamaterials with a simple exper- 
imental setup that does not require additional filtering, 
amplifiers or biasing of the sample. This experimental 
simplicity together with the robustness and tunability of 
the evanescent signal over a wide range of incident pow- 
ers and temperatures make these experimental results all 
that more compelling. Most noteworthy, however, is the 
unanticipated transition from one set of resonances to a 
second set of resonances. This transition could be a man- 
ifestation of the arrays transitioning from being locked in 
phase to being unlocked; however, further investigations 
are needed to confirm this speculation or provide an al- 
ternative explanation. 
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